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' The fully hadronic final states of two signal processes from an mSUGRA inspired sce- 

nario (SUSY-P5) are studied within a full simulation of the LDC' detector model. These 
are chargino pair and neutralino pair production, i.e. e + e~ — > xtxt ~* 1<j'Xi 99' Xi 
and e + e~ — > X2X2 ~~ > 99Xi 99Xi- Both processes have to be separated sufficiently from 
all background to measure the respective production cross sections and extract the 
masses of the involved bosinos, m(xt), m{5&) an d friLSP = m (xl)- This is achieved by 
fitting the energy spectra of the reconstructed gauge bosons while taking into account 
the finite width of the boson mass. From simulation data corresponding to 500 fb _1 of 
r*^ luminosity, a mass resolution of about 0.5 GeV seems to be achievable. 

<D 

^! 1 Introduction 

^ , 

Due to the clean environment and precisely known inital state at the ILC, it will be possible 
to measure very small signal cross sections and extract the masses of SUSY particles even 
from fully hadronic final states. As an example, from the SUSY-P5 scenario [2], pair pro- 
duction of charginos (xi Xi) an( i neutralinos (X2X2) an d their subsequent hadronic decays 
I/*-) I (<?<7Xi 99~Xi) are studied. To successfully extract the masses of intermediate SUSY particles 

0^ ' from a fully hadronic final state, it is mandatory to reduce the background to a minimum. 

Especially difficult is the mutual separation of both signal processes. Four jets need to be 
assigned to two gauge bosons, either two W or two Z bosons, depending on whether they 
originated from a decaying chargino or neutralino. In addition, a considerable amount of 
missing energy and momentum carried away by two of the lightest supersymmetric particles 
(LSPs) impede a complete reconstruction of the final state. 

All Monte Carlo samples are part of the DESY mass production [3] . The events are gener- 
ated with Whizard [4], run through the full simulation of the LDCPrime_02Sc detector model 
and finally reconstructed with MarlinReco [5] using Pandora Particle Flow (PFlow) [6]. 



2 Selection strategies and kinematic fits 

This study concentrates on the reduction of SM background using a few very basic require- 
ments and a kinematic fit to improve separation of the chargino from the neutralino signal. 
A rough background scan as depicted in Fig. 1 (a) shows WW and WWZ production to 
be the dominant source of SM background. Especially, 6-fermion final states from WWZ 
production, with an invisible Z-decay, i.e. qq'qq'vD, constitute a nearly irreducible back- 
ground component. However, the processes shown in Fig. 1 are not normalised to the same 
luminosity. While the SUSY signal processes correspond to 500 fb _1 , background samples 
were only available for « 20-50 fb^ 1 for QCD multijet production (qq, qqqq final states) and 
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1 fb 1 for 77/67 processes. Background from other SUSY processes has also not yet been 
considered, but it is expected to not render the preliminary results of this study obsolete. 




/ \ 1 000 2000 3000 4000 5000 /i \ °° 70 80 90 100 110 120 130 140 

(a) dst file # | (b) W energy [GeV] 



Figure 1: (a) SM background contributions. Events left after requiring: black: Ntracks > 20 
and 150 GeV < E ms < 300 GeV; blue: lepton veto with Ei ept on < 20 GeV; red: E jet > 5 GeV 
and | cos(0jet)| < 0.9 for all jets, (b) Empirical fit to the energy spectrum of the W boson 
consisting of a 3 rd -order polynomial convoluted with a Gaussian (see text). 

Some very basic requirements reduce leptonic events and events from multijct QCD pro- 
duction (see caption of Fig. 1), while the very similar neutralino signal is hardly affected. To 
separate the chargino and neutralino signals from each other, two hypotheses are compared 
for the invariant di-jet masses, which can either form a W or a Z boson. The resulting 
X 2 -values for each hypothesis are plotted against each other, so that the final selection cri- 
terium for chargino events results in: 

X 2 (W, Jlj2 ) + X 2 (W, hh) < 2 with X 2 (W, j k ji) = (m w - m(j k , ji)f /(5 GeV) 2 
X 2 (Z, Jlj2 ) + X 2 {Z, j 3 j 4 ) > 4 with X 2 (Z, j k jt) = (m z - m( 3klJl )) 2 /(5 GeV) 2 . 

A more detailed description of the entire procedure is given in [7] . After the above selection 
of chargino-type events, the energy spectrum of di-jet masses in Fig. 1 (b) largely resembles 
that of the W boson a . The spectrum is fitted with an empirical function consisting of a 
3 rd -order polynomial (4+2 par. to describe the edge positions) convoluted with a Gaussian 
(2 par.) yielding the following values for the bosino masses: 

m(xf) = 219.7 ± 1.54 GeV <-> 216.5 GeV (nominal chargino mass) 
m(x?) = 121.2 ± 0.59 GeV <-> 115.7 GeV (nominal ncutralino/LSP mass). 

The central fit values differ rather much from the nominal bosino masses in the SUSY-P5 
scenario. However, for an analysis of real data, such a shift can be compensated by tuning 
the MC distribution such that the experimental data is well reproduced. 

A different approach postpones the reduction of events from SM background in favour of 
a slightly better separation of chargino and neutralino signals with the help of a 2-dim. cut 
in the (VI , V2) mass plane, with V being a W or Z boson (see Fig. 2). The cut is tuned, but 
not yet optimised, for a high efficiency and purity of the chargino pair production (e x tt) . 

a SM background is multiplied by 10, since so far the luminosity is only about 50 fb -1 for most background 
processes. It is also not included in the fit, while the neutralino pair background is taken into account. 



LCWS/ILC2008 



A simple kinematic lC-fit, requiring equal di-jet masses (EMC), is applied to the three 
di-jet combinations of all events that can be clustered into 4 jets. The fit is meant to help 
find the correct jet-pairings, improve the mass and energy resolutions for the gauge bosons 
and thus, ultimately, the starting point for a fit to the W (Z) boson energy spectra. The 
di-jet masses populating the (VI, V2) mass plane in Fig. 2 correspond to the jet-pairings 
with the best fit probability and are assumed to be the correct ones. 
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Figure 2: Selection of chargino-type events: all events between the middle (red) and the 
lower (orange) lines are assumed to originate from chargino decays, while events between 
the middle and upper lines are taken to be neutralino decays. 



3 Physics-driven fit to the energy spectrum of the W boson 



Having selected the events between the lower two lines in Fig. 2 as chargino events, the 
energy spectrum of the corresponding gauge boson is then fitted with the masses of the 
SUSY particles as free parameters. However, this time, not an empirical fit function is 
used but a physics-driven function. Firstly, the gauge boson's energy and momentum are 
expressed in terms of the three relevant masses mi = m(xx, X2), mi — m(xx) — m LSP, 
and my = raw z- 
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From this the inverse expression has to be calculated for cither the chargino (neutralino) 
mass. Next, the effect of the finite boson mass width is added by convolusion with a Breit- 
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is added to the emerging fit function: 
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Oj : jet energy resolution . 
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Figure 3: Fit to the energy spectrum of the W boson (a) on generator level and (b) after 
reconstruction and a first rough chargino selection. 



The fit values for the reconstructed bosino masses are actually quite close to the nominal 
generator level values: 

m(xf) — 215.3 ± 2.45 GeV «-> 216.5 GeV (nominal chargino mass) 
TO (x?) = H5.2 ± 0.46 GeV <-> 115.7 GeV (nominal ncutralino/LSP mass). 

Although the error on the chargino mass is still large, this is expected to be reduced by 
a more sophisticated selection (higher e x ir) and a simultaneous fit to the cross sections of 
both signal processes, i.e. chargino and neutralino pair production. The cross sections could 
then serve as a further input to the fit of the gauge boson energy spectra. 

While only the dominant SM background (WWZ — > qq'qq'vv) is considered for the 
distributions in Fig. 3 (a,b), signal and background processes are normalised to the same lu- 
minosity of 500 fb _1 . This provides an estimate of the sensitivity of the chargino (neutralino) 
selection to a possible contamination from SM background. 



4 Conclusions and outlook 

Physics analyses using a full detector simulation play an important role in the ongoing 
optimisation effort for the ILC detector. This study, although performed on simulation data 
of a previous detector model (LDCPrime_02Sc), constitutes one part of the general effort for 
the ILD detector concept. In cases where the gauge boson masses are reconstructed from 4- 
jet final states, the mass and energy resolution is expected to improve drastically due to the 
interplay between a highly segmented calorimeter with tracking ability for neutral particles 
and the new approach of Particle Flow algorithms. Together with more sophisticated fits to 
the resulting gauge boson energy spectra, mass resolutions of 0.5-1 GeV are achievable for 
those SUSY particles accessible to the ILC (xf, X°j an d X?)- 

Next steps will be to include all SM and relevant SUSY background, optimise the selection 
and fit procedures and utilise different polarisation scenarios to improve the S/B-ratio. 
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